ACEEE Int. J. on Signal & Image Processing, Vol. 01, No. 03, Dec 2010 



An Optimized Transform for ECG Signal 

Compression 

Mrs.S.O.Rajankar 1 , Dr. S.N. Talbar 2 

'Sinhgad College of Engineering, Pune, India. 

Email: supriya.rajankar@gmail.com 

2 S.G.G. S. College of Engineering and Technology, Nanded, India. 

Email: sntalbar@yahoo.com 



Abstract-^ significant feature of the coming digital era is the 
exponential increase in digital data, obtained from various 
signals specially the biomedical signals such as 
electrocardiogram (ECG), electroencephalogram (EEG), 
electromyogram (EMG) etc. How to transmit or store these 
signals efficiently becomes the most important issue. A digital 
compression technique is often used to solve this problem. 
This paper proposed a comparative study of transform based 
approach for ECG signal compression. Adaptive threshold is 
used on the transformed coefficients. The algorithm is tested 
for 10 different records from MIT-BIH arrhythmia database 
and obtained percentage root mean difference as around 
0.528 to 0.584% for compression ratio of 18.963:1 to 23.011:1 
for DWT. Among DFT, DCT and DWT techniques, DWT has 
been proven to be very efficient for ECG signal coding. 
Further improvement in the CR is possible by efficient 
entropy coding. 

Index Terms: Electrocardiogram (ECG), Discrete Fourier 
Transform (DFT), Discrete Cosines Transform (DCT), 
Discrete Wavelet Transform (DWT) etc. 



I. 



Introduction 



The ECG is a one of the important physiological signal 
which depicts the electrical activity of a heart. ECG 
processing is a topic of great interest in the scientific 
community because based on the ECG's a diagnosis is 
done for detecting abnormalities in the heart functioning 
[7]. The amount of ECG data grows with the increase of 
number of channels, sampling rate, sample resolution, 
recording time, etc. As an example, with a sampling rate of 
360 Hz, 11 bits/sample data resolution, a 24-h record 
requires about 43 Mbytes per channel[5],[6]. Therefore an 
effective data compression scheme for ECG signal is 
required in many practical applications such as ECG data 
storage, ambulatory recording systems and ECG data 
transmission over telephone line or digital 
telecommunication network for telemedicine. The main 
goal of electrocardiogram (ECG) data compression 
techniques is to preserve the most useful diagnostic 
information while compressing a signal to an acceptable 
size. Lossless compression is the best choice as long as the 
compression ratio is acceptable, but it cannot usually offer 
a satisfactory compression ratio (CR). To obtain significant 
signal compression, lossy compression is preferable to a 
lossless compression. 

ECG data compression algorithms can be categorized as 
direct data compression (DDC), Transform based and 
Parameter extraction. Direct data compression is a time 
domain compression algorithm which directly analyses 
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samples where inter-beat and, intra-beat correlation is 
exploited. Some of the DDC algorithms are Turning point 
(TP), Amplitude zone time epoch coding (AZTEC), Co- 
ordinate reduction time encoding system (CORTES), Scan 
along polynomial approximation (SAPA), and Non 
redundant template. These algorithms suffer from 
sensitiveness to sampling rate, quantization levels and high 
frequency interference. It fails to achieve high data rate 
along with preservation of clinical information [10]. In 
Transform based technique compressions are accomplished 
by applying an invertible orthogonal transform to the 
signal. Due to its de-correlation and energy compaction 
properties the transform based methods achieve better 
compression ratios. The commonly used transforms for 
data compression are DCT, DST, DFT and DWT. In 
parameter extraction methods a set of model 
parameters/features are extracted from the original signal 
(model based) which involves methods like Linear term 
prediction LTP) and analysis by synthesis (ASEC) [5], 
[10]. 

In this paper different transform like DCT, DFT and 
DWT are studied for ECG signal compression. A threshold 
based algorithm is proposed to achieve better compression. 

II. Transform Techniques 

A. Discrete Fourier Transform 

Discrete Fourier Transform is a fundamental transform 
in digital signal processing with applications in frequency 
analysis, signal processing etc. The periodicity and 
symmetry properties of DFT are useful for compression. 
The w* DFT coefficient of length N sequence [f(x)} is 
defined as in "(1)": 



F{u) = Y^f{x)e- n ' ruxlN u = 0,l,....N -1 



And its inverse transform (IDFT) as in "(2)": 



(I) 



/(*) = —$>("> 



jlituxlN 



N 



0,1, ,N-1 



(2) 



The number of complex multiplications and additions to 
compute DFT is N~. Moreover fast algorithms exist that 
makes it possible to compute DFT efficiently [9]. This 
algorithm is popularly known as Fast Fourier Transform 
(FFT) which reduces the computations to N log 2 N. 
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B. Discrete Cosine Transform 

Discrete Cosine Transform is a basis for many signal 
and image compression algorithms due to its high de- 
correlation and energy compaction property. A discrete 
Cosine Transform of A' sample is defined as in"(3)": 



F(u) 



2 K r 7i(2x + V,u 1 

C(u)2 J f(x)cos[ — :r7r - 1 



2N 



u = 0,l,...,N-l, 

Where C(u) 



(3) 



for u=0 



1, otherwise. 



The function fix) represents the value of xth samples of 
input signals. F(u) represents a DCT coefficients. The 
inverse DCT is defined in similar fashion as follows in 

"(4)": 



ti ^ Pv'rr m?< s r x(2x + l) U l 

/(■*) = J — / C(u)F(u)cos[ ] 

\JN f^ 2N 

x = 0A,...,N -I. 



(4) 



Since DCT belongs to family of DFT there are Fast 
DCT algorithms of computational complexity N log 2 N 
similar to FFT[9]. 

C. Discrete Wavelet Transform 

Wavelet transforms have become an attractive and 
efficient tool in many applications especially in coding and 
compression of signals because of multi-resolution and 
high-energy compaction properties. Wavelets allow both 
time and frequency analysis of signals simultaneously 
because of the fact that energy of wavelet is concentrated in 
time and still possesses the wave like characteristics. As a 
result wavelet representation provides a versatile 
mathematical tool to analyze transients, non-stationary 

signals. Wavelets y/ a b (t) are functions generated by one 
single basis function, called mother wavelet y/ (t) by 
dilation 'a' and translation 'V represented as in"(5)": 



Va, b tt) 



1 



w 



t-b 



v 



a 



(5) 



If a and b are discredited then discrete wavelets can be 
represented as in"(6)": 



V / j, k ( t ) = a o 1 ¥{a~'t-kb ) 
where for dyadic sampling a =2 and £> =1 

y/ jk {t) = 2- ill y/{2- i t-k) j,k^Z 



(6) 



(7) 



©2010 ACEEE 
DOL01.IJSIP.01.03.224 



34 



This multi-resolution wavelet algorithm decomposes a 
signal /(7)with the help of (/){t ) and wavelet function 
by y/ (t). These functions together resolves the signal into its 

coarse and detail components[8],[10]. Thus using multi- 
resolution idea, the signal f(t)\& defined in terms of scale 
and wavelet coefficients, c(k) and d(k), respectively. That 



is. 



f(t) = £ c(*)&(0 + £ £ d{j,k)ys. k (t) 



fc=-oo 



j=Q fe=-oo 



(8) 



The first summation gives a function that is a low 
resolution or a coarse approximation of f(t); the second one 
represents the higher or finer resolution to give detailed 
information of the signal. 

III. PERFORMANCE MEASUREMENT 

The evaluation of performance for testing ECG 
compression algorithms includes three components: 
compression efficiency, reconstruction error and 
computational complexity. The compression efficiency is 
given by compression ratio (CR). The compression ratio 
and the reconstruction error are usually dependent on each 
other. The computational complexity component is part of 
the practical implementation consideration [6], [8]. 

A. Compression Measurement 

All data compression algorithms minimizes data storage 
by reducing the redundancy wherever possible, thereby 
increasing the compression ratio. The compression ratio 
(CR) is defined as the ratio of the number of bits 
representing the original signal to the number of bits 
required to store the compressed signal. A high 
compression ratio is typically desired. A data compression 
algorithm must also represent the data with acceptable 
fidelity while achieving high CR. 



CR = 



original 



compressed 



(9) 



B. Distortion Measurement 

One of the most difficult problems in ECG compression 
and reconstruction is defining the error criterion that 
measures the ability of the reconstructed signal to preserve 
the relevant information. Different objective error measures 
namely; root mean square error (RMSE), percentage root 
mean difference (PRD), signal to noise ratio (SNR) are 
used for calculation of reconstruction error. Among these 
error measures, the ECG processing is frequently measured 
by the percentage root mean difference (PRD). It is most 
commonly defined as: 



PRD 



X„=i x <» -*(«) 



E>« ! 



xlOO 



(10) 
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where x(n) and x(n) are the values of the original and 
reconstructed samples, respectively, and N is the length of 
the window over which the PRD is calculated. The PRD 
indicates reconstruction fidelity by point wise comparison 
with the original data. Despite their wide use, PRD do not 
indicate precisely the quality of signal's reconstruction and 
the decompressed signal has to be evaluated by visual 
inspection also. 

IV. CODING ALGORITHM 

Transform based techniques because of their high 
compression ability have gained popularity. In this method 
the preprocessed signal is transformed to get the de- 
correlated coefficients. The thresholding or quantization of 
transformed coefficients gives the actual compression, 
which is lossy one. But it has good performance and low 
computational cost. In this paper the comparative study of 
DFT, DCT and DWT is carried out. Following are some 
important steps of the proposed algorithm. 

• Transform the ECG signal using DFT/DCT/DWT. 

• To achieve an adaptive threshold compute the 
maximum value of the transformed coefficients. 

• Apply the threshold of a fix percentage based on 
absolute maximum values of the transform 
coefficients. 

• After threshold operation to reconstruct the signal 
apply inverse transform. 

For wavelet based compression the signal is 
decomposed up to fifth level using Bior4.4 wavelet. Bi- 
orthogonal wavelets also allow perfect reconstruction of 
the signal using linear phase filter banks which in turn 
avoid reconstruction error at the beginning and at the end 
of the signal. 

V. EXPERIMENTAL RESULTS 

The data from MIT-BIH arrhythmia database is used to 
test the performance of the coding schemes. This database 
is sampled at 360Hz, and the resolution of each sample is 
11 bits/samples. The algorithm is tested for 1024, 2048 & 
4096 samples of the record numbers 100, 103, 105, 107, 
111, 200, 203, 205, 207 and 210. Table I summarizes the 
results for record number 207 with maximum 
threshold. "Fig. 1" shows original and reconstructed 



waveform of the ECG signal. The CR and PRD values 
increases proportionally if the percentage threshold is 
increased as in "fig. 2" and "fig.3".The energy compaction 
ability for different transform is shown in "fig.4". With 
4096 samples, for DFT method the PRD is around 0.222 
to 1.020% with the corresponding CR 4.362:1 to 
15.814:l.For DCT the PRD is around 0.181 to 0.9346% 
with CR 5.446:1 to 19.23:1. In case of DWT the PRD is 
around 0.214 to 0.584% with the CR 9.315:1 to 23.011:1. 
The "fig. 5" explains about effect of number of samples on 
CR. 
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Figure 1. (a) Original and (b) Reconstructed signal of record number 207 
with DWT (PRD=0.584% and CR=23.01 1:1) 
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Figure 2. Comparison of CR for different transforms 
(record number 207). 



TABLE I. 
PERFORMANCE RESULTS FOR DIFFERENT TRANSFORM METHODS FOR RECORD NUMBER 207 WITH MAXIMUM THRESHOLD 



Number of 
samples 


DFT 


DCT 


DWT 


CR 


PRD 


CR 


PRD 


CR 


PRD 


1024 


10.343 


0.405 


12.337 


0.172 


18.963 


0.528 


2048 


11.702 


0.683 


14.027 


0.623 


21.558 


0.561 


4096 


15.815 


1.021 


19.230 


0.935 


23.011 


0.584 
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Figure 3. Comparison of PRD for different transforms 
( record number 207). 
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Figure 4. Comparison of % energy retained for different transforms 
(record number 207) 
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Figure 5. Comparison of CR for different number of samples with DWT 
(record number 207) 



Conclusion 

For the transforms considered in this paper as the 
percentage threshold is increased the number of 
transformed coefficients with zero value gets increased 
which improves the CR but PRD also gets increased 
proportionally. As the PRD indicates reconstruction fidelity 
the increase in its value is actually undesirable. In 
comparison with DFT and DCT the rise in PRD for DWT 
is less but CR is more, also energy compaction in DWT 
more. So the DWT is found to be more suitable for ECG 
signals compression. Also it is observed that as the number 
of samples is increased the correlation among the samples 
and hence the CR is increased. It is possible to further 
increase the CR by using effective coding methods. 
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